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Using updated gluon distributions from global fits to data, we investigate the sensitivity of direct 
photoproduction of heavy quarks and exclusive production of vector mesons to varying strength of 
gluon modifications. Implications of using these processes for constraining nuclear gluon distribu- 
tions are discussed. 
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Ultraperipheral relativistic heavy ion collisions can ex- 
plore several aspects of particle and nuclear physics and 
have been extensively discussed in the literature (for a 
small sample of references, see e.g. flT-[l2j]). In this arti- 
cle, we investigate the sensitivity of direct photoproduc- 
tion of heavy quarks and exclusive production of vector 
mesons to different gluon distribution functions. This 
idea, originally proposed in Ref. Q , can be used to con- 
strain the possible distribution functions from data on 
production of heavy quarks and of vector mesons. Here 
we report on the study of certain distribution functions, 
not considered so far. 

A key ingredient of our calculations, the photon flux in 
ultraperipheral collisions, can be evaluated by the equiv- 
alent photon (Weizsacker- Williams) method [l]. Im- 
provements of the method has been documented in sev- 
eral publications For a gi ven impact parameter 
b, the flux of virtual photons with photon energy k 
is d 3 N^(k,b) /dkcPb, depending strongly on the Lorentz 
factor 7. At the Large Hadron Collider (LHC) at CERN 
the Lorentz jl factor in the laboratory frame is 7455 for 
p-p and 2930 for Pb-Pb collisions. The photon flux also 
depends strongly on the adiabaticity parameter £ = kb/j 
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Due to the modified Bessel functions Kq(() and Ki(Q, 
the photon flux possesses the asymptotic property of 
an exponential drop-off at £ > 1, above a cutoff en- 
ergy determined essentially by the size of the nucleus, 
Ecutoff ~ 7 MeV jb (fm). The relationship between 
the Lorentz contraction factor associated with the rela- 
tive velocity between the colliding nuclei, and the col- 
lider energy per nucleon, E/A, in GeV, is given by 
7 = 2 7 £ - 1 « 2(1.0735£/yl) 2 . 

Integrating d 3 N~ f (k,h)/dkd 2 b over impact parameters 
with the constraint of no hadronic interactions yields the 
total photon flux dNj(k)/dk. An analytic expression for 
this total flux, strongly dependent on the reduced adia- 
baticity parameter — 2kR,A/j for AA collisions, or 
(ft A = k(R p + Ra)M f° r pA collisions, is derivable in the 
approximation whereby Eq. ([1]) is integrated over impact 
parameters larger than the sum of the radii of the partic- 
ipants. While this is a good approximation, a relatively 



better estimate of the total flux is obtained by taking the 
average over the target surface d, Hi"! 
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Although both the analytic expression and Eq. ([2]) , evalu- 
ated numerically, have been utilized in the present study, 
we report only results using the numerical flux. The dif- 
ferences between the results from analytic and numerical 
fluxes are generally of the order of 10 — 15%. With the 
knowledge of the photon flux, any generic total photopro- 
duction cross section can be factorized into the product 
of a photonuclear cross section cr x (k) and the photon 
flux, dNj/dk, a x = J dk(dN 7 /dk)a x (k). 

The photonuclear processes described in the present 
work are dependent on gluon distributions in nuclei. This 
dependence influences the structural characteristics of 
these processes, especially in the case of the exclusive 
photoproduction of vector mesons, in which the gluon 
distribution enters quadratically. It is a rather well- 
known fact that the distributions of partons (i.e. quarks 
and gluons) in nuclei are quite different from the distri- 
butions in free nucleons, due to the complex, many-body 
effects in the nuclear medium. This is expected to mani- 
fest in experimental observables such as the cross section 
and rapidity distributions. 

Nuclear parton distribution functions, F^(r,x,Q 2 ), 
are often for technical convenience expressed as a con- 
volution of "nuclear modifications" Ra( r , x, Q 2 ) and free 
nucleon parton distribution functions f a (x, Q 2 )- Here the 
subscript a denotes a parton species and the superscript 
A a particular nucleus. The variables are the position 
vector r, parton momentum fraction x (Bjorken-x), and 
a hard scale (factorization scale) Q 2 . However, since lim- 
ited availability of data does not permit a determination 
of the spatial dependence, current nuclear parton distri- 
butions from global fits are functions of x and Q 2 only. 
The nuclear effects encoded in the nuclear modifications 
Ra(x, Q 2 ) can be categorized based on different intervals 
in x. At small values of x (x < 0.04), we have the phe- 
nomenon generally referred to as shadowing. This is a 
depletion, in the sense that in this interval, the nuclear 
parton distributions are smaller compared to the corre- 
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sponding distributions in free nucleons, i.e. < 1. 

Antishadowing, which is an enhancement (R^ > I), oc- 
curs in the range 0.04 <jr < 0.3. Another depletion, 
the classic EMC effect [13|], is present in the interval 
0-3 ^ x < 0.8, while for x > 0.8, the Fermi motion region, 
we have another enhancement. It is important to note 
that although both shadowing and the EMC effect (an- 
tishadowing and Fermi motion) correspond to depletion 
(enhancement), the physical principles and mechanisms 
governing these phenomena are quite different. Further 
details can be found in 

The determination of gluon distributions in both nu- 
cleons and nuclei is a highly nontrivial task. In the 
usual determination of parton distributions from global 
fits to data, a preponderance of the experimental data 
is from Deeply Inelastic Scattering (DIS) and Drell-Yan 
(DY)processes. Since gluons are electrically neutral, 
their distributions cannot be directly extracted from DIS; 
they are inferred from sum rules and the Q 2 evolution of 
sea quarks distributions. The situation is even worse in 
the nuclear case: the available data is much less than for 
nucleons, and there is the added complication of a mass 
dependence. It is therefore not unusual for nuclear gluon 
distributions from different global fits to differ signifi- 
cantly, especially in the magnitude of the various nuclear 
effects (shadowing, antishadowing, etc). This is espe- 
cially obvious at low Q 2 (i.e. around their initial starting 
scales) since evolution to high Q 2 tends to lessen the dif- 
ferences. Earlier global analyses [f8l - [2ll ] relied heavily on 
fixed-target nuclear deep-inelastic scattering (DIS) and 
Drell-Yan (DY) lepton-pair production, with the atten- 
dant low precision and weak constraints on nuclear gluon 
distributions. The use of extended data set, incorporat- 
ing data on inclusive hadron production in deuteron-gold 
collisions, has been pioneered in 22, 23], with better con- 
straints on gluon modifications. We should also mention 
the approach in [24| which utilizes the Gribov picture of 
shadowing. Despite all these advances the nuclear gluon 
distribution is still the least constrained aspect of global 
fits to nuclear parton distributions, as large uncertainties 
still persist at both small and large x. 

Four recent gluon distributions are utilized in the 
present study. For the nucleon gluon distributions we 
use the Martin-Stirling-Thorne- Watts (MSTW08) par- 
ton distributions (25|. In the nuclear case we use three 
nuclear modification sets. Two sets are by Eskola, 
Paukunnen, and Salgado, namely EPS08 and EPS09 [H, 
HI]. The third is the Hirai-Kumano-Nagai (HKN07) dis- 
tributions [U. The gluon distributions from MSTW08 
serve two purposes: as the free nucleon distributions used 
in conjunction with nuclear modifications, and also as a 
" special" nuclear gluon distribution in the absence of nu- 
clear effects. The latter usage is particularly useful for 
highlighting the influence of the various nuclear effects on 
observables. Our calculations are to leading order (LO); 
thus all distributions are evaluated at this order. 

These four gluon distributions have different character- 
istics due to the different strength of their nuclear modi- 



fications. As already stated, one can view the MSTW08 
gluon distributions as nuclear gluon distributions in the 
limit of zero nuclear effects ( R^ = I). In Fig. [I] we 
show the nuclear modifications from EPS08, EPS09, and 
HKN07 at the factorization scale Q 2 = M 2 ^, appropri- 
ate for the elastic photoproduction of the J/vf meson. 
At this scale, HKN07 has a rather weak gluon shadow- 
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FIG. 1: (Color Online) Nuclear gluon modifications in Pb, 
Rg b {x,Q 2 = Mj /9 ), from EPS08 (dashed line), EPS09 (solid 
line), and HKN07 (dash-dotted line) respectively. 

ing which extends well into the antishadowing region, no 
antishadowing and gluon EMC effect, and an early onset 
of Fermi motion. EPS09 exhibits a moderately strong 
shadowing, and appreciable antishadowing and EMC ef- 
fect, with a quite strong Fermi motion. Nuclear modifica- 
tions are strongest in EPS08: an especially strong shad- 
owing, and substantial antishadowing, EMC, and Fermi 
motion. Thus in terms of shadowing we have a progres- 
sion from zero effects to weak effects, moderate (inter- 
mediate) effects, then to strong effects as one progresses 
from MSTW08 to EPS08. We have not included uncer- 
tainties from gluon distributions in this study. Further 
discussions on uncertainties in nuclear parton distribu- 
tions can be found in (23j. 

From the point of view of a Fock space decomposi- 
tion, photon interactions with hadrons and nuclei can 
be direct or resolved. In direct interactions the photon 
behave as a pointlike particle ("bare photon") while in 
resolved interactions the photon fluctuates into a quark- 
antiquark state or an even more complex partonic state. 
In the present study we will focus attention only on the 
direct photon contribution. The contribution of resolved 
photon processes in ultrapcriphcral heavy ion collisions 
is treated in detail in Ref. [8j. 
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Photon-gluon fusion leading to the production of a 
heavy quark pair is the dominant subprocess when high 
energy photons are incident on a nucleus. Due to the 
high energies involved, perturbative QCD is applicable, 
and the cross section can be expressed as a convolution 
of the partonic cross section for the subprocess — > qq 
and the nuclear gluon distribution: 



^A^qqX 



(W lA ) 



dx er' 1 



\W ig )G A {x,Q 2 m 



with x the momentum fraction carried by the gluon, and 
x m in = 4TOq/W^ A . Here m q is the mass of the heavy 
quark (charm or bottom), W lA (Wig) denotes the center- 
of-mass energy of the photon- nucleus (photon-gluon) sys- 
tem, and G A (x,Q 2 ) is the nuclear gluon distribution. 
The pQCD factorization scale, Q, is quite arbitrary; the 
cross section is thus scale dependent. In the present study 
we use two different scales in order to assess the magni- 
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tude of scale dependency: a dynamic scale Q — . , , ; 
as in (2^, and a static scale Q 2 — Am 2 (m 2 ) for charm 
(bottom) respectively, as in [8j. We take m c = 1.4 GeV 
and m b =4.75 GeV for consistency with the MSTW08 
parton distributions. 

The cross section for the photon-gluon fusion leading 
to quark pair production is [26l428j : 



2na em a s {Q 2 )e 2 



W 2 



(1 + p - \p 2 ) ln(^=) - (1 + 



with e q the electric charge of the quark, a em the elec- 
tromagnetic coupling constant, and /? = 4m 2 /W 2 g . The 
strong coupling constant, a s (Q 2 ), needed for the calcu- 
lation, is evaluated to one loop at the scale Q 2 using the 
evolution code contained in the MSTW08 package. The 
total photoproduction cross section ct(A[7]A — > AqqX) 
is obtained by convoluting the equivalent photon flux, 
dN 7 (k)/dk with a ^ A ^ x (k): 



a(A[ 7 ]A -> AqqX) = I dk a^ A ^ x (k) 
' dk 



(5) 



The final state qq rapidity is dependent on the pho- 
ton energy k and the gluon momentum fraction x: x = 
(W lg /W lA )e v . Thus changing variable from k to y and 
differentiating (da/dy = kda/dk), the differential cross 
section with respect to rapidity is thus da lA ^"^ x j dy = 
[kdN 1 {k)ldk)a^ A ^ qliX {k). 

Both the total cross section and the rapidity distri- 
bution involve the product of the photon flux and the 
photonuclear cross section. The flux decreases exponen- 
tially with increasing photon energy k while the cross 
section, due to its dependence on x m i n , increases with 
k since x m i n is inversely related to k. This interplay 
not only decisively influences the magnitude of the to- 
tal cross section and rapidity distributions, but also the 



relative contributions of the various x-dependent nuclear 
modifications of the gluon distribution in G A (x, Q 2 ). For 
instance, gluon shadowing is of less importance in photo- 
production of bottom quarks compared to charm quarks. 
This is due to the the fact that the advent of shadowing 
contribution in 66 production occurs at a larger photon 
energy where the flux is more suppressed than in cc pho- 
toproduction where the onset occurs at a lower photon 
energy, and thus with a larger usable photon flux avail- 
ability. 

The differential cross section for the elastic photo- 
production of a vector meson V on a nucleus A in 
the exclusive process A [7] A —> AAV can be writ- 
ten as da-> A ^ VA /dt = da~> A -> VA /dt\ t=Q \F(t)\ 2 , where 

da~* A ^ VA /dt\ t =o is the forward scattering amplitude and 
F(t) is the nuclear form factor. The dynamical informa- 
tion is encoded in the forward scattering amplitude while 
the momentum transfer of the elastic scattering is deter- 
mined by the form factor, which is dependent on the 
spatial attributes of the target nucleus. 

Diverse mechanisms have been employed in the eval- 
uation of the dynamical content of the forward scatter- 
ing amplitude for heavy mesons. In this study we use 
the simple amplitude calculated from leading order two- 
gluon exchange in perturbative QCD [29|, |3fJ and cor- 
rected for other relevant effects (such as relativistic cor- 
rections, inclusion of the real part of the scattering am- 
plitude, next-to-leading order NLO effects, etc, see for 
instance [HIHH) through a phenomenological multiplica- 
tive correction factor (y . For elastic photoproduction on 
protons, the corrected LO scattering amplitude can be 
written as: 



dt 



t=0 



16ir 3 a 2 s T e 
3aM v 



■)[xg p (x,Q 2 )] 2 . (6) 



Here, My is the mass of the vector meson (J/vl/ and 
T(ls) in the present study), x = M 2 /W 2 p is the frac- 
tion of the nucleon momentum carried by the gluons, and 
g p {x, Q 2 ) is the gluon distribution in a proton, evaluated 
at a momentum transfer Q 2 — (My/2) 2 . Eq. [6] is easily 
generalized to the nuclear case: 



da^ A ~ 



■VA 



dt 



Cv 



16n 3 a 2 s T e 
3aM y 



[xG A (x,Q 2 )} 2 , (7) 



where G A (x,Q 2 ) = g p (x,Q 2 
gluon distribution and R A (x, Q 



xR A (x,Q 2 ) 



is the nuclear 
the gluon modification. 
The correction factor £y is estimated by constraining 
the calculated cross sections for elastic vector mesons 
photoproduction on protons, a n,p ^ yVp (W 7P ), to reason- 
ably reproduce the photoproduction data from HERA: 
[111 for J/# and (aJUl for T(ls). a ~< p ^ Vp is obtained 
through a ^ Vp = (1/6) da^ Vp /dt\ f=() with slope pa- 
rameter 6. Using 6 = 4.5 GeV -2 , we have = 1/3.5 
and Or (is) = 1-0, and the results are displayed in Fig. [2] 
Or (is) = 1-0 for the elastic photoproduction of T un- 
derestimates older HERA data [33, w4| but seems to be 
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FIG. 2: (Color Online) Cross section for photoproduction of 
J/ty (T) as a function of energy W lp . In the upper panel 
the solid line is the LO result while the dashed line is the 
corrected LO result using = 1/3.5, with data is from 

[33| . In the lower panel data taken from [33t |35| and solid 
line depicts LO result. 



FIG. 3: (Color Online) Forward scattering amplitude for elas- 
tic photoproduction of J/SP in PbPb collisions as a function 
of momentum fraction x for four different gluon distributions. 
Dotted line depicts result using the MSTW08 gluon distri- 
bution (no nuclear modifications). Dashed, dot-dashed, and 
solid lines are results from nuclear-modified gluon distribu- 
tions from EPS09, EPS08, and HKN07 parton distributions 
respectively. 



adequate for the newer analysis in (35[. This could indi- 
cate the smallness of higher-order effects and other cor- 
rections, whereas they are of major importance in J/^T/ 
production. 

In contrast to photoproduction of heavy quarks, the 
quadratic dependence of the differential cross section on 
the gluon distribution has the significant implication of 
making exclusive vector meson production a very sensi- 
tive probe of nuclear gluon modifications. This is appar- 
ent from Fig. [3] where the forward scattering amplitude 
for «//*& production in PbPb collisions at LHC energy 
is plotted as a function of x for the four gluon distribu- 
tions under consideration. The different characteristics 
displayed in Fig. Q] are clearly manifested. 

The nuclear form factor, F(t) 1 is given by the Fourier 
transform of the nuclear density distribution: F (t) = 
J d 3 r p{r) e lq r , where q is the momentum transferred. 
For a heavy nucleus it is customary to model the density 
distribution as a Woods-Saxon distribution with param- 
eters from electron scattering, p(r) = po/[l + e^ r 
with central density p$, radius Ra and skin depth d. For 
208 Pb in use at the LHC, p = 0.16/fm 3 , R A = 1.2A 1 ' 3 
fm, and d — 0.549 fm [36j. Since the Fourier transform 
of a Woods-Saxon density distribution does not admit of 
an analytic form, we employed the commonly-used mod- 
ified hard sphere (a convolution of a hard sphere with a 



Yukawa term) 0,0, H3] to approximate p(r) in F(t) 
4-irpo 



F{q 



\t\) 



Aq 3 



sm(qRA) — qRA cos(qRA) 



1 



a 2 q 2 



(8) 



The range of the Yukawa term, a, is 0.7 fm, and the form 
factor is a simple product of the Fourier transforms of 
the hard sphere and the Yukawa term. 

The photonuclear cross section is thus given by 



(T 



da^ A ~ 



■VA 



dt 



t=0 Jt mi „(k) 



dt\F(t)f (9) 



Here t m j n (k) = (M 2 /4/c7l) 2 , as is appropriate for narrow 
resonances [40j | . The total cross section is a convolution 
of the photonuclear cross section and the photon flux: 



A[~i]A^AAV _ 



dk dN^l ulA ^VA {k) 
dk 



dk 



dN 7 (k) 
dk 



dt- 



da"i A - 



■VA 



dt 



\F(t)\ 2 . (10) 



t=o 



It is often of practical interest to represent the cross 
section in terms of the rapidity of the vector meson. 
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The photon energy, fc, is related to the rapidity, y, by 
k = (My/2) exp(y). Using this relationship, the differ- 
ential cross section with respect to rapidity is given by 
da~< A ^ VA /dy = {kdN 1 (k)/dk)cr' lA ^ VA {k).' Thus with a 
knowledge of the photon flux the differential cross sec- 
tion, da/dy, is a direct measure of the vector meson pho- 
toproduction cross section for a given photon energy. 

We now discuss the results of our calculations for both 
the inclusive photoproduction of heavy quarks (cc and 
bb) and the exclusive production of vector mesons (J/^f 
and T(ls)). For the cross sections for heavy quarks we 
present results for the two scales employed: a dynamic 



scale Q — s and a static scale Q = Am 2 c (Q 2 



777 ; 



;) for 



charm (bottom), while for the vector mesons we present 
both the LO and corrected LO results. In the case of 
rapidity distributions we show only the result for the dy- 
namic scale for heavy quarks, and the corrected LO for 
the vector mesons. Since we deal with symmetric PbPb 
collisions, both nuclei can serve as source/target and the 
total rapidity distribution is the sum of both, and sym- 
metric about y = 0. 
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FIG. 4: (Color Online) Rapidity distributions of cc photopro- 
duction in PbPb collisions at the LHC for photons incident 
from the right (upper panel) and from the left (lower panel). 
Solid line depicts result using the MSTW08 gluon distribu- 
tion (no nuclear modifications). Dashed, dot-dashed, and 
dotted lines are results from nuclear-modified gluon distribu- 
tions from EPS09, EPS08, and HKN07 parton distributions 
respectively. 

The sensitivity of heavy quark photoproduction to nu- 
clear gluon modifications is more transparent in rapid- 
ity distributions than in total cross sections. In Fig. 0] 
we show the rapidity distributions for cc production in 
ultraperipheral PbPb collisions at the LHC, employing 



the four gluon distributions described earlier. The upper 
panel depicts the distributions with the incident photon 
coming from the right, while the lower panel shows the 
distributions with the incident photons from the left. The 
total, which is the sum of both panels, is displayed in 
Fig. El 
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FIG. 5: (Color Online) Total rapidity distributions of the pho- 
toproduction of cc in PbPb collisions at the LHC. Solid line 
depicts result using the MSTW08 gluon distribution (no nu- 
clear modifications). Dashed, dot-dashed, and dotted lines 
are results from nuclear-modified gluon distributions from 
EPS09, EPS08, and HKN07 parton distributions respectively. 

As mentioned earlier, the value of x m in(k) determines 
the nuclear effects contributing to the photonuclear cross 
section a lA ^ qqX (k). Since the rapidity distribution 
da/dy is directly proportional to this cross section, it 
is rather straightforward to elucidate the features of the 
rapidity distribution based on the relative contributions 
of the relevant nuclear effects. Fig. Q] is thus quite helpful 
in understanding the characteristics of the distributions 
displayed in Fig. 0] and, by extension, Fig. [5] 

Let us start from midrapidity (y = 0) and move to- 
wards higher positive values of y. In the upper panel of 
Fig. @] (photons from right) the contribution from shad- 
owing becomes progressively smaller, such that at around 
y = 3 shadowing is negligible and antishadowing becomes 
the dominant contributor. At around y — 6 antishad- 
owing peters out, and the distributions are subject only 
to EMC/Fermi motion effects. On the other hand, for 
the lower panel in Fig. @] (photons from left), shadowing 
becomes progressively stronger and more dominant for 
increasing y. The effect of shadowing in this region is 
mitigated by the rather strong suppression of the photon 
flux at large fc; the distributions thus rapidly die out with 
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TABLE I: Total cross sections for direct photoproduction of 
cc in ultraperipheral PbPb collisions at the LHC. All cross 
sections are in millibarns (mb). 



Gluon Distribution 


Q* = s 
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MSTW08 


1170 


1090 


EPS08 


890 


780 


EPS09 


1000 


910 


HKN07 


1080 


1000 



increasing positive y. The opposite trend is observed for 
increasingly negative y values starting from midrapidity. 
Here, for photons incident from the right (upper panel), 
increasing negative values of y implies increasing shad- 
owing and flux suppression. For photons incident from 
the left (lower panel), the transition is from shadowing 
to EMC/Fermi motion effects. 

We now consider the characteristics of the total ra- 
pidity distributions due to the amalgamation of the con- 
trasting tendencies exhibited in the upper and lower pan- 
els of Fig. |U Shadowing is the dominant nuclear effect 
for —3 < y < 3, and the rapidity distributions in this 
region reproduce the observed trend of gluon shadow- 
ing strength as exhibited in Fig. [Q MSTW08 with its 
zero gluon shadowing gives the largest rapidity distri- 
bution while EPS08, with its strong gluon shadowing, 
gives the smallest. Due to strong flux suppression, shad- 
owing is most markedly apparent for the rapidity range 
— 2 < y < 2. This range therefore provides a good win- 
dow to discriminate among different gluon shadowing sce- 
narios. 

The rapidity intervals 3 < y < 6 corresponds to x m i n 
in the antishadowing region (deep shadowing) for right 
(left) incident photons and vice versa for —6 < y < —3. 
Due to the photon flux suppression in the deep shadow- 
ing region, the rapidity distributions are sensitive mainly 
to antishadowing in addition to both EMC effect and 
Fermi motion. Since both EPS08 and EPS09 have sub- 
stantial antishadowing, their rapidity distributions reflect 
this, being slightly larger than those from HKN07 and 
MSTW08. The discriminatory power here is not as ap- 
preciable as in the shadowing case though, due largely to 
the smallness of the distributions. 

For both rapidity ranges y < —6 and y > 6, x m i n > 0.2 
and the relevant contributing nuclear effects are the EMC 
and Fermi motion since the contribution from antishad- 
owing is small, and that from shadowing practically 
nonexistent by virtue of flux suppression. Due to the be- 
havior of HKN07 in this interval (no EMC effect, only en- 
hancement), da/dy from HKN07 is largest. Both EPS08 
and EPS09 nuclear modifications exhibit EMC effect and 
Fermi motion, and the destructive interference from both 
effects render their rapidity distributions to practically 
coincide with that from MSTW08. 

In Table [J we present the total cross section for the di- 
rect photoproduction of cc at two different scales as dis- 
cussed earlier. For both scales the cross sections exhibit 
a clear trend in conformity with the relative strength of 



gluon shadowing. MSTW08 with no modifications gives 
the largest cross section while EPS08 gives the smallest, 
due to its strong gluon shadowing. It is also noteworthy 
that the static scale gives lower cross sections relative to 
the dynamic scale. The difference between the two scales 
increases progressively with increasing shadowing, from 
fa 7% for MSTW08 to around 13% for EPS08. Overall, 
the total cross section seems also a good discriminator of 
different gluon shadowing scenarios. 

Our cross sections can be compared with results from 
other studies on cc photoproduction. In [lj| [3j| the cross 
section with no shadowing (equivalent to our MSTW08) 
is 1250 mb, from EKS98 [ijfy (somewhat akin to EPS09) 
it is 1050 mb, while from FGS (strong shadowing as in 
EPS08) [II it is 850 mb. Likewise in U the no shad- 
owing cross section is 1790 mb while the cross section 
from EKS98 is 1500 mb. The study in [9| gives the cross 
section as 2056 mb. Overall, our results are closest to 
those reported in [l(| [39[ . Differences in results are not 
only attributable to the different gluon distributions and 
photon fluxes used in the earlier studies, but also pos- 
sibly to different values of the running strong coupling, 
a s{Q 2 ), which enters multiplicatively in the expression 
for the photon-gluon cross section. 

We now discuss our result for the total cross sections 
and rapidity distributions for bb production. The cross 
sections are orders of magnitude less than in cc produc- 
tion, and also correspondingly exhibit less sensitivity to 
nuclear modifications. In Fig. [6] we show the rapidity 
distribution for bb in ultraperipheral PbPb collisions at 
the LHC. Shadowing dominates in the rapidity interval 
— 2 < y < 2 and is most clearly manifested in the ra- 
pidity window — 1 < y < 1. Thus this interval presents 
the best sensitivity to shadowing effects in bb produc- 
tion. Although less marked, the progression of rapidity 
distribution with relative shadowing strength follows the 
trend observed in cc production: MSTW08 still gives the 
largest distribution while EPS08 gives the smallest. 

As in the case of cc production, there is a slight mani- 
festation of the influence of antishadowing around y = —3 
and y = 3. The distributions practically overlap for 
y < —4 and y > 4; thus overall, the interval — 1 < y < 1 
seems to afford the best sensitivity to nuclear effects, in 
this case primarily shadowing. More detailed treatment 
of the x-dependence of bb production at the LHC is pre- 
sented in [381 ]. 

TABLE II: Total cross sections for direct photoproduction of 
bb in ultraperipheral PbPb collisions at the LHC. All cross 
sections are in millibarns (mb). 



Gluon Distribution 


Q A = i 


Q 1 =mt 


MSTW08 


6.2 


7.0 


EPS08 


5.8 


6.2 


EPS09 


6.0 


6.6 


HKN07 


6.1 


6.7 



The total photoproduction cross sections for bb from 
the four gluon distributions are presented in Table |IIJ It 
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FIG. 6: (Color Online) Total rapidity distributions of bb pho- 
toproduction in PbPb collisions at the LHC. Solid line de- 
picts result using the MSTW08 gluon distribution (no nuclear 
modifications). Dashed, dot-dashed, and dotted lines are re- 
sults from nuclear-modified gluon distributions from EPS09, 
EPS08, and HKN07 parton distributions respectively. 



is readily observed from the table that nuclear modifi- 
cations have manifestly weaker effects on bb production 
relative to cc, since the results from two widely different 
distributions like MSTW08 and EPS08 are almost in the 
same ballpark. This observation is also apparently scale 
independent. Unlike the cc case though, the static scale 
gives larger cross sections than the dynamic scale, and the 
difference decreases with relative shadowing strength. 

Our cross sections can be compared with results from 
other studies on bb photoproduction. In [HI H^j the cross 
section with no shadowing is 4.9 mb, from EKS98 it is 
4.7 mb, while from FGS it is 4.4 mb. Likewise in [1] the 
no shadowing cross section is 0.718 mb while the cross 
section from EKS98 is 0.686 mb. The study in [j| gives 
the cross section as 20.1 mb. Again, our results are closest 
to the values reported in [Iol |. The comment concerning 
the differencies in the cc results is also applicable here. 
Additional details can be found in 0. 

It is pertinent at this point to remark on the limita- 
tions inherent in our calculated cross sections for both 
cc and bb photoproduction. The results are to leading 
order; thus higher-order effects have not been taken into 
account either explicitly or through a phenomenological 
correction factor. In addition, we have not included the 
resolved contributions which are quite sizeable (see Q). 
We thus advocate that these limitations should be borne 
in mind, moreso in view of the disparities in cross section 
results from the present work and previous studies. 



We now present our results on elastic photoproduction 
of the J/vP and T(ls) in the framework of a leading- 
order two-gluon exchange formalism in QCD. Apart from 
the quadratic dependence on the gluon distribution, two 
other quantities namely the integrated nuclear form fac- 
tor and the photon flux affect the attributes of both the 
rapidity distribution and total cross section. The photon 
flux has support at low photon energy k, which translates 
to large x, while the integrated form factor favors large 
k, or equivalently, small x. 

Fig. [7] shows the rapidity distribution of the J/'F in 
ultraperipheral PbPb collisions at the LHC. Shadowing 
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FIG. 7: (Color Online) Total rapidity distributions of exclu- 
sive photoproduction of J in PbPb collisions at the LHC 
in the modified hard sphere density distribution approxima- 
tion. Solid line depicts result using the MSTW08 gluon distri- 
bution (no nuclear modifications). Dashed, dot-dashed, and 
dotted lines are results from nuclear-modified gluon distribu- 
tions from EPS09, EPS08, and HKN07 parton distributions 
respectively. 

is the relevant nuclear effect in the rapidity interval 
— 3 < y < 3 and unsurprisingly, the rapidity distributions 
mimic the behavior in the shadowing region of Fig. [T] 
The largest rapidity distribution is given by MSTW08, 
followed by HKN07, and EPS09. The smallest is by 
EPS08 due to its strong gluon shadowing. The rapid- 
ity window — 2 < y < 2 manifestly depicts the signifi- 
cant distinction between the various gluon distributions 
arising from the quadratic dependence. Antishadowing 
manifests in the intervals — 5 < y < — 4 and 4 < y < 5; 
the effect though is quite slight. 

Table IIIII shows the total cross sections for the elastic 
photoproduction of the J /\& using the four gluon distri- 
butions considered in our study. The leading order result 
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TABLE III: Total cross sections (in mb) for elastic photopro- 
duction of J/ty in ultraperipheral PbPb collisions at the LHC. 
Second column is the result from LO calculation while third 
column displays result from corrected LO calculation. 



Gluon Distribution 


LO 


Corrected LO 


MSTW08 


260 


74 


EPS08 


36 


10 


EPS09 


101 


29 


HKN07 


173 


49 



is under the second column while the corrected leading 
order result is under the third. As expected, the total 
cross sections reproduce the trend seen in the rapidity 
distributions: MSTW08 gives the largest cross section 
while EPS08 yields the smallest. 

There have been studies of the photoproduction of J/^P 
in ultra-peripheral collisions at LHC using diverse pro- 
duction mechanisms 0, 0> EH 0, EHil]. Here, since 
we are interested in the sensitivity to gluon modifica- 
tions, we use the simple leading order (LO) two-gluon 
exchange mechanism corrected for additional relevant ef- 
fects through a multiplicative factor. Our work is similar 
in spirit to the study reported in [l2j in the sense that 
four different gluon distributions were also utilized. Thus 
we compare our results to the work in [l2j which reports 
a no shadowing cross section of 74 mb, a cross section 
from EPS08 of 13 mb, EKS98 of 39 mb, and from DS03 
[13] (somewhat similar to HKN07) of 61 mb. While the 
corrected LO results are somewhat close to these values, 
the uncorrected LO results presented in Table IIIII are 
consistently higher. 

In Table IIVI we also present the equivalent cross sec- 
tions for J /"J/ production in ultraperipheral AuAu colli- 
sions at RHIC. As is readily apparent, the results follow 
the shadowing trend as observed for PbPb collisions at 
the LHC. 

TABLE IV: Total cross sections (in fib) for elastic photopro- 
duction of J/ty in ultraperipheral AuAu collisions at RHIC. 



Gluon Distribution 


LO 


Corrected LO 


MSTW08 


1222 


349 


EPS08 


699 


200 


EPS09 


868 


248 


HKN07 


902 


258 



The rapidity distributions for T(ls) production is 
shown in Fig. [FJ and exhibit identical trend observed for 
J/ty. Shadowing remains the relevant nuclear modifi- 
cation for practically the entire rapidity range shown in 
the figure, and is markedly manifested in the interval 
— 2 < y < 2. Thus rapidity distribution in this inter- 
val should be a good discriminator of gluon shadowing 
strength. 

In Table [V] we show the total cross sections for the 
elastic photoproduction of the T(ls) for the four gluon 
distributions under study. Unsurprisingly, the total cross 
sections reflect the progressive trend of the relative shad- 
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FIG. 8: (Color Online) Total rapidity distributions of exclu- 
sive photoproduction of T(ls) in PbPb collisions at the LHC 
in the modified hard sphere density distribution approxima- 
tion. Solid line depicts result using the MSTW08 gluon distri- 
bution (no nuclear modifications). Dashed, dot-dashed, and 
dotted lines are results from nuclear-modified gluon distribu- 
tions from EPS09, EPS08, and HKN07 parton distributions 
respectively. 



TABLE V: Total cross sections for elastic photoproduction of 
T(ls) in ultraperipheral PbPb collisions at the LHC. 

Gluon Distribution Cross Section (fib) 

MSTW08 189 

EPS08 99 

EPS09 130 

HKN07 146 



owing strength, MSTW08 giving the largest cross section 
while EPS08 yields the smallest. 

Previous studies of the photoproduction of T(ls) in 
ultra-peripheral collisions at LHC energies have been re- 
ported in [H HI El IH|. As in the case of J/* pro- 
duction, we compare our results with the values reported 
in p~2t | . The no shadowing cross section is 163 fib, from 
EPS08 22 fib, EKS98 120 ^b, and from DS03 148 fib. 
Except for the case of EPS08, our results presented in 
Table fvl are quite close to these values. 

In conclusion, we have considered the direct photopro- 
duction of heavy quarks (charm and bottom) and elastic 
photoproduction of vector mesons (J/'S and T(ls)) in 
ultraperipheral PbPb collisions at LHC energy. These 
two processes are dependent on nuclear gluon distribu- 
tions, and are therefore potentially useful in constraining 
modifications such as shadowing and antishadowing in 
nuclear gluon distributions. In order to assess the sensi- 
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tivity to these modifications, we have utilized four recent 
gluon distributions, chosen on the basis of the relative 
strength of their modifications. For each process we con- 
sidered two observables: rapidity distributions and total 
cross sections. 

In direct photoproduction of heavy quarks the gluon 
dependence is linear and different modifications are su- 
perimposed due to the integration over the momentum 
fraction x. Despite these, rapidity distributions for cc 
manifest appreciable sensitivity to shadowing around 
midrapidity and a slight sensitivity to antishadowing at 
more forward and backward rapidities. Thus cc photo- 
production offers good constraining potential for shad- 
owing, and a somewhat less potential for antishadowing. 
Although photoproduction of bb is less sensitive to mod- 
ifications than cc, the influence of shadowing is evident 
around midrapidity, and it thus offers some constraining 
ability for shadowing. 

The quadratic dependence on gluon modifications 



makes elastic photoproduction of vector mesons particu- 
larly attractive for constraining purposes. This is man- 
ifestly apparent from the rapidity distributions for both 
J and T(ls) photoproduction which exhibit very good 
sensitivity to gluon shadowing over an appreciable range 
about midrapidity. Thus both offer remarkable poten- 
tial in constraining the shadowing component of nuclear 
gluon distributions. 

Determination of nuclear modifications from ab-initio 
calculations of cross sections is beset with difficulties. A 
more feasible approach is to compare photoproduction 
in proton- nucleus and nucleus-nucleus collisions, where 
many theoretical uncertainties and systematic errors can- 
cel (see [13] )■ Further work along this line is in progress. 

We acknowledge support by the US Department of 
Energy Grants DE-FG02-08ER41533 and DE- FC02- 
07ER41457 (UNEDF, SciDAC-2) and the Research Cor- 
poration. 
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